Phonon-assisted current noise in molecular junctions 
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We investigate the effects of phonon scattering on the electronic current noise through nanojunc- 
tions using the non-equilibrium Green's functions formalism extended to include the counting Held. 
In the case of weak electron-phonon couphng and a single broad electronic level we derive an analytic 
expression for the current noise at arbitrary temperature and identify physically distinct contribu- 
tions based on their voltage dependence. We apply our theory to the experimentally relevant case 
of a D2 molecule placed in a break-junction and predict a significant inelastic contribution to the 
current noise. 
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Introduction. — The fabrication of atomically sharp con- 
tacts has opened up the possibility of creating junc- 
tions formed by a single molecule bridging metal elec- 
trodes However, the inherent complexity of this 
emerging field poses fundamental challenges and makes 
junctions formed by very simple molecules (e.g. hydro- 
gen 0, S, i] or water [|) an invaluable testbed from 
both the experimental and the theoretical point of view. 
In these systems, inelastic effects due to the interaction 
between transport electrons and molecular vibrational 
modes (phonons) result in an abrupt change of the differ- 
ential conductance at the onset of the phonon emission. 
These features have been exploited to establish unam- 
biguously the presence of the molecule in the contact 13 
and, when combined with shot-noise measurements 0,0], 
they allow for a detailed characterization of the junction. 

Theoretical descriptions of inelastic transport through 
nanojunctions have so far focused mainly on the current- 
voltage characteristics, while less attention has been paid 
to the study of noise. Phonon-scattcring effects on the 
differential conductance have been addressed both with 
i^S ll|Jl^ and with simphfied 



ab-initio methods [7 
(one-level) models 



16l |. Noise calculations 
based on one-level models have also been put forward 
within the rate equation approach [isl . [l^ or within the 
non-equilibrium Green's functions (NGF) formalism [l3| 
with a mcan-ficld-like approximation for the noise. In 
this work we study inelastic effects on the current noise 
with the NGF approach taking consistently into account 
all the correlations due to phonon-assisted scattering up 
to a given order in the e-ph interaction. We apply our 
theory to the case of a hydrogen-bridge junction and 
predict a significant inelastic contribution to the current 
noise. 

Model & methods. — The system we consider can be 
represented as a central device region which is tunnel- 
coupled to non-interacting metallic leads H = He + 
Hl.r + Ht- Neglecting the spin degree of freedom, 
the central region can be described by the Hamiltonian 



He = Ho + J2e fi^ib% + Et Mi^dldj{bl+be), where 
d] and arc the electron and phonon creation operators, 
Ho = Hq is the single-particle effective Hamil- 
tonian of the electrons moving in a static arrangement 
of atomic nuclei, and is the e-ph coupling matrix for 
the £-th phonon mode. Here, boldface notation stands 
for matrix over electronic space. The leads and tunneling 
Hamiltonians are given by Hl,r = Y.k.a=L,R ^a,kC^a,k^a.k 

and Ht = Efc,a=L,i?(^a,fccl,fcrfi -f h.c). The states in 
the leads are occupied according to the Fermi distribu- 
tion /„(e) = [1 + e/Jf^-A"-)]-!^ where [3 = l/keT is the 
inverse temperature and fia is lead-a chemical potential, 
with — Mi? = fixed by the applied bias voltage. 

To calculate the current and the noise we employ the 
extended NGF technique [3| to find the cumulant gen- 
erating function S{X) = 1/ (ze)A -I- S/{2e^)X^ + • • • , from 
which the current / and the noise 5" can be calculated 
straightforwardly. The key idea [l3| is to modify the 
Hamiltonian by adding a time-dependent phase A(t)/2 
to the tunneling matrix elements VI with X{t) = ±A 
for t on the upper and the lower Keldysh contour, respec- 
tively. It has been shown in Ref. [20| for the Anderson 
model, that d\S{X) is related to the Keldysh Green's 
function Gx{t,t') = ~i{Tcd{t)d\t'))x, where the expec- 
tation value is now evaluated in the presence of A(t). 
Generalizing that result to the multilevel case, we write 



TT{TL[e-^\l - fL)G-+ + e^'fLG+-]} 



dSjX) 
dX J 2TTh 

(1) 

Here rL^ (e) = 2^Efe " ^fc,") is the level 

broadening due to the coupling to lead-o: and G^"*", 
are appropriate components of the solution of the 
non-equilibrium Dyson equation G\ = gA + gA^AGA. 
The check sign indicates matrices in the Keldysh 
space and the superscripts —/+ correspond to the for- 
ward/backward branch of the Keldysh-contour. The ma- 
trix gA is Green's function of the system in the presence 
of the leads and of the counting field but without the 
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e-ph interaction. Its inverse is given by 



el- 
-iTLe 



- Ho - i Ea Ta [fa (s) - 1/2] zFlC^ (s) + iTRfR{e) 

-'^[1 - hie)] - iTr[1 - fnie)] -e 1 + Ho - z ^M<.{e) - 1/2] 



(2) 
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As gA already includes the coupling to the leads, is 
the self-energy solely due to the e-ph coupling. Being in- 
terested in the weak coupling limit, we expand the Dyson 
equation to the lowest (second) order in the e-ph coupling 



gA + Ex^^x^gx, where T:)^' is given by the Fock 



(2) 



diagram (77, rj = ±) 



de' 



(£)-*E / ^df(e-e')M,gf(£')M,. (3) 



The Hartree term has been neglected since it cannot 
contribute by any truly dynamical features in which we 
are primarily interested Above, cQ^{e) stand for free 
phonon Green's functions for the £-th phonon mode 
df^ie) = E.=± [ - Mm + me + Shu:,) ± T'^] 

and dj'^{e) ^ -2m[{Mt + l)5{e ±niui>) + Nt5{eT^t)\, 
with Ml the (generally non-equilibrium) occupation of 
mode 

Truncating the Dyson equation to the second order 
in directly yields the expressions for Ga=o and 
9AGA|^_g, which are the ingredients to evaluate the cur- 
rent and the noise. Integration over energy can be per- 
formed analytically assuming the electronic structure to 
be slowly changing over few multiples of a typical phonon 
energy around the Fermi level Ep and approximating 
r„(e) « T^{Ef) and g^^o(e) « e\=q{Ef) 1, i E E | . 

Physically important non-equilibrium phonon heating 
effects [m can be qualitatively taken into account by a 
rate equation for the average phonon occupation num- 
ber Ml 0, which can be viewed as a kinetic-equation- 
like approximation to the full NGF studies [13, [3- 
In the broad-level approximation introduced above, this 
leads to a bias-dependent occupation number Mi[V) = 
nsi^i) + adniiy), where nsiuj) = l/{e^'^'^ — 1) is the 
Bose distribution, ad is a parameter which depends on 
strength of the external phonon damping, and ni{V) 
takes into account the power dissipated by the trans- 
port electrons into the phonon mode Q. In the fol- 
lowing we will focus on the two opposite regimes of (i) 
thermally equilibrated phonons {ad = 0), and (ii) non- 
equilibrated phonons [ad = 1). In this case it is niiV) ~ 
{\eV\/huJi - l)e{\eV\ - huji)/A for keT < huji 

Analysis. — For sake of clarity, we focus here only on 
the case of a single electronic level eo with symmetric cou- 
pling to both leads Vl ~ F/j = F and coupled to a single 
phonon mode with frequency woj occupation Mq and cou- 
pling constant M . Already such a simple model reveals 



many essential features of inelastic transport through 
nano-j unctions such as the phonon-induced step behav- 
ior of the differential conductance 0, [IB] • Before focus- 
ing on the inelastic corrections to the noise, we shortly 
reexamine the results for the current. We find that the 



current through the device is given by / = I^i - 
/ci = {e^/h) X TV and 



1, with 



27r 



(1 - 2T) 



W{V -l)-W{V + 1) 



+ (2AAo + l)(3-4r)y , 



(4) 



with the reduced voltage V = eV/hujQ, the dimen- 
sionless e-ph couphng 7eph = M^T^/F^, and W{x) = 
xcoih-{l3huiQx/2). Here T is the elastic transmission co- 
efficient r = 1G"^12f2 = FV(A2 + f2) and A = {Ep-eo) 
gives the position of the single level with respect to the 
Fermi energy. The inelastic current results then from the 
sum of two contributions with distinct behavior with re- 
spect to the bias voltage: while the first term of Eq. ([1]), 
which is responsible for the step features in the non-linear 
conductance, saturates to constant values for \V\ > 1, the 
second one grows linearly with V for Mo = nsitOQ) and 
quadratically in the case of non-equilibrated phonons. 
This second term has a clear physical interpretation 
in terms of electrons experiencing the coupling to the 
phonon as a stochastic quasi-static shift in the energy 
of the level. This in turn affects the transmission coeffi- 
cient, which becomes dependent on the displacement of 
the oscillator T{Q). Averaging over Q and retaining only 
terms to the second order in M one obtains 



(TiQ))-^ 



(A-MQ)2 + r2 



■^«T + 7eph(3-4r)(g2), 



(5) 

where (•) indicates the average over a stationary distri- 
bution of the oscillator, so that (Q) = and (Q^) = 
{2M0 + 1). We can therefore interpret the last term of 
Eq. ([4]) in terms of an elastic-like contribution with av- 
eraged transmission over the fluctuating position of the 
oscillator [13] 

We now turn to the main result of our work which 
is the phonon-assisted noise. The current noise through 
the device is given by S' = Sd + Si^id, with the standard 
expression for the elastic noise S",,! = {e"^ /h){T'^ / f3+T{l — 
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FIG. 1: (Color online) Inelastic noise Sinci as a function of the 
bias voltage at zero temperature T = for different values of 
the transmission coefficient. Both cases of equilibrated {aa = 
0, thick lines) and non-equilibrated phonons (od = 1, dashed 
lines) are shown. Inset: Plot of ci |{r^o,|eV|=fit^o} (^^U line) 
and (ci + «;i)|{T=o,ey=o} (dash-dotted line) as a function of 
T (A/o = for these parameter values). 
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sW'iV + s) 



s=±l 



(6) 



being the correction due to inelastic scattering. Here 
W'{x) = dW/dx and the coefficients read cq ~ 4(2A/'o + 
1)T(5 - er), ci = (2AAo + l)(12r2 - UT + 3), Kq = 
4T(1 - r)(l + 2A/-o){ [1 + 2ns(^o)] - 2/f}hu;o} - (1 - 2T)^ 
Ki = 4r(l - r)[l + 2nB(c^o)] - (1 - 2r)2(l + 2AAo), 
and, finally. K2 = 2T(1 — 2T). In the zero bias limit 
Sinoi satisfies the fluctuation-dissipation theorem, i.e. 
S'inci|y=o = 2(?inci//3, where C7i„ci = (9/i„ci/W)|y=o is 
the inelastic correction to the linear conductance. 

Similarly as for the current, we interpret the first term 
of Eq. ([6]) as a quasi-static correction to the elastic noise 
due to averaged transmission over the displacement of 
the oscillator. In fact, along the same line which lead to 
Eq. (I5|), it is easy to see that 7ephCo and 7ephCi correspond 
exactly to the contribution of order to {T{Q)'^) and 
{T{Q){1 — T{Q))), respectively. The remaining terms 
are dynamic contributions which take into account the 
phonon exchange effects. Analogously to the current, the 
quasi-static contribution has a distinct voltage behavior 
compared to the dynamic one for large voltage, being the 
only one which does not saturate for \V\ > 1. 

In the limit of zero temperature Eq. ^ simplifies no- 
ticeably becoming 



5i 



inol T=0 



ki|t=o • (1 



{ko|t=o + ci|t=o • 1^1 

\v\)eii~\v\)}. 



(7) 



FIG. 2: (Color online) Temperature dependence of inelastic 
noise, (a) Sinei for the case of equilibrated phonons (a^ = 
0) at fcflT = O.OSfttJo (a typical experimental value 0, Hi)- 
Different lines correspond to various values of the transmission 
coefficient T (color code identical to Fig. [T|. (b) Same as 
in (a) but at higher temperature ksT = 0.25?ia;o. Inset: 
Plots of h{y\ = kl^iV + 1) - 2W{V) + W{V - 1)] and 
f2{V) = W'{V + 1) - W'{V - 1) at kgT = 0.03 hujo (left 
panel) and ksT = 0.25/icJo (right panel). These functions 
characterize the dynamic contributions to Sind in Eq. (|6]). 



The inelastic noise is then a piecewise function charac- 
terized by the coefficient ci|t=o for \V\ > 1 and by the 
combination (ci + ki)\t=o for \V\ < 1, which takes into 
account the competition between quasi-static contribu- 
tion and the dynamic one. The signs of these coefficients 
determine whether Sinci is an increasing or decreasing 
function of voltage, leading to the rich behavior presented 
in Fig.[TJ Interestingly, depending on the value of T both 
(ci -I- ki)|t=o and ci|t=o can be negative, resulting in 
a negative contribution to the inelastic noise. In other 
words, for a rather wide range of transmissions, phonon 
scattering events lead to suppression of the current noise 
through the device. Note that both ci and ki may de- 
pend on voltage via Aq. However, since the phonon heat- 
ing becomes effective only above the phonon emission 
threshold, S'inci is always a linear function for \V\ < 1, 
while energy accumulation into the phonon mode results 
in the quadratic increase of S'inci for \ V\ > 1 in the case of 
no external damping (non-equilibrated phonons, ad = 1). 

For finite temperatures qualitatively new features ap- 
pear (see Fig. [5]). For small temperatures the curves get 
rounded around V = 1 and for high transmission even an 
additional dip occurs (see Fig. [2^). The changes become 
more pronounced if the temperature is of the order of 
the phonon frequency, when the kinks of S'inci are largely 
washed out. Nevertheless, it is still possible for a wide 
set of parameters to have a negative inelastic correction 
to the noise and even a sign change at some finite \V\ 
(see Fig. [35). 

Results. — We apply now our formulas to the case of a 
single hydrogen molecule between platinum atomic con- 
tacts [ij. Experimentally, it has been shown that hydro- 
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FIG. 3: (Color online) Total current noise S — Sd + ^inci 
through a D2 molecule and its derivative dS/dV as a func- 
tion of voltage. The elastic contributions Sd and dSd/dV 
axe plotted as dashed lines for comparison. Parameter values 
fcsT = 0.029 fttJo, htoo = 50meV, r = 0.9825, jeph = 0.011 
and ad — 1 are taken from Ref. and correspond to the 
conductance measurements in Ref. The noise level at 

eV = hwo = 50 meV corresponds to experimentally accessi- 
ble 3 X 10^^'^A'^/Hz. Phonon heating {ad ~ 1) responsible 
for the finite slope of dS/dV for \eV\ > tkuo must be included 
because of the large mass mismatch between the D2 molecule 
and Pt atoms. 



gen can form a stable bridge between Pt electrodes with 
conductance close to the quantum unit carried dom- 
inantly by a single, almost transparent channel Such 
a picture has been confirmed by first principle calcula- 
tions showing that a single conductance channel forms 
due to strong hybridization between the H2 anti-bonding 
state and the Pt metal states, while the bonding state is 
not involved in the transport [23| . 

Fig. [3] represents our prediction for the phonon-assistcd 
noise through a D2 junction, where typical values for T, 
ksT and 7eph have been taken from Ref. Q and cor- 
respond to the experimental data of Ref. 0] . One main 
message of our work is that, despite of the very weak e-ph 
coupling, inelastic corrections give a sizable contribution 
to the total noise through a D2 junction for \eV\ > tkjQ. 
On the other hand, inelastic corrections are negligible 
for \eV\ < hwQ, thus justifying the interpretation of noise 
measurement in this regime in terms of elastic theory (as 
done in 

In conclusion, we have presented a perturbative scheme 
for the calculation of the inelastic contribution to the cur- 
rent noise in systems with weak e-ph interaction. In the 
experimentally relevant case of a single broad level, we 
have derived an analytic expression for the inelastic noise 
at arbitrary temperature and distinguished terms that 
correspond to simple renormalization of the transmission 
coefficient from those which contain true dynamical ef- 
fects. Applying our theory to the case of a D2 junc- 
tion, we predict a sizable contribution to the total noise 
due to inelastic processes. Our scheme can be straight- 
forwardly extended beyond the present model to cases 
with multiple electronic levels and phonon modes, asym- 



metric coupling to leads, energy-dependent transmission, 
and/or moderate e-ph coupling with application in cur- 
rent ab-initio methods 0, [3, [13, 11 1 . 
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